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Abstract 

It is shown that the classical L-operator algebra of the elliptic Ruijsenaars-Schneider 
model can be realized as a subalgebra of the algebra of functions on the cotangent bundle 
over the centrally extended current group in two dimensions. It is governed by two 
dynamical r and f-matrices satisfying a closed system of equations. The corresponding 
quantum R and i?-matrices are found as solutions to quantum analogs of these equations. 
We present the quantum L-operator algebra and show that the system of equations on 
R and R arises as the compatibility condition for this algebra. It turns out that the 
i?-matrix is twist-equivalent to the Felder elliptic .R^-matrix with R playing the role of 
the twist. The simplest representation of the quantum L-operator algebra corresponding 
to the elliptic Ruijsenaars-Schneider model is obtained. The connection of the quantum 
L-operator algebra to the fundamental relation RLL = LLR with Belavin's elliptic R 
matrix is established. As a byproduct of our construction, we find a new TV-parameter 
elliptic solution to the classical Yang-Baxter equation. 



*Stcklov Mathematical Institute, Vavilov 42, GSP-1, 117966, Moscow, Russia; arut@class.mi.ras.ru 
T Steklov Mathematical Institute, Vavilov 42, GSP-1, 117966, Moscow, Russia; chekhov@class.mi.ras.ru 
*Stcklov Mathematical Institute, Vavilov 42, GSP-1, 117966, Moscow, Russia; frolov@class.mi.ras.ru 



1 



1 Introduction 



The appearance of classical dynamical r-matrices |T], || in the theory of integrable many-body 
systems raises an interesting problem of their quantization. On this way one may hope to 
separate the variables explicitly. 

At present, the classical dynamical r-matrices are known for the rational, trigonometric 0, [3] 
and elliptic [|], || Calogero-Moser (CM) systems, as well as for their relativistic generalizations 
- the rational, trigonometric || [7[] and elliptic || || Ruijsenaars-Schneider (RS) systems [|TO|] . 
It is recognized that dynamical systems of the Calogero type can be naturally understood in the 
framework of the Hamiltonian reduction procedure O, |T2|. Moreover, the reduction procedure 
provides an effective scheme to compute the corresponding dynamical r-matrices [[3], |l3f . 

Depending explicitly on the phase variables, the dynamical r-matrices do not satisfy a single 
closed equation of the Yang-Baxter type, that makes the problem of their quantization rather 
nontrivial. In [f[4}| , the spin generalization of the Calogero-Sutherland system was quantized 



by using the particular solution [15] of the Gervais-Neveu-Felder equation |16, 15] and in 



[T7|1 , it was interpreted in terms of quasi- Hopf algebras. This system is not integrable, but zero- 



weight representations of the quantum L-operator algebra admit a proper number of commuting 
integrals of motion. However, it seems to be important to find such a quantum L-operator 
algebra for the Calogero-type systems that possesses a sufficiently large abelian subalgebra. 

Recently, an algebraic scheme for quantizing the rational RS model in the i?-matrix formal- 
ism was proposed [n|. We introduced a special parameterization of the cotangent bundle over 
GL(N,C). In new variables, the standard symplectic structure was described by a classical 
(Frobenius) r-matrix and by a new dynamical f-matrix. The classical L-operator was intro- 
duced as a special matrix function on the cotangent bundle. The Poisson algebra of L inherited 
from the cotangent bundle coincided with the L-operator algebra of the rational RS model. It 
is this reason why we called L the L-operator. Quantizing the Poisson structure for L, we found 
the quantum L-operator algebra and constructed its particular representation corresponding to 
the rational RS system. This quantum algebra has a remarkable property, namely, it possesses 
a family of N mutually commuting operators directly on the algebra level. 

It is well-known that the elliptic RS model is the most general among the systems of the 
CM and RS types. In this paper, we are aimed to include this model in our scheme. Recall 
that the classical L-operator algebra for the elliptic RS model can be obtained by means of 



the Poissonian ||T9[ or the Hamiltonian reduction schemes POfl . In the first scheme, the affine 
Heisenberg double is used as the initial phase space and in the second one, the cotangent 
bundle over the centrally extended current group in two dimensions is considered. Thus, the 
appropriate phase space we choose to deal with the elliptic RS model is the cotangent bundle 
T*GL(N)(z, z) over the centrally extended group GL(N)(z, z) of double loops. The application 
of our approach |18|] is not straightforward since one should work with the infinite-dimensional 
phase space and, therefore, the correct description of the Poisson structure on T*GL(N)(z, z) 
in the desired parameterization requires an intermediate regularization. 

Describe briefly the content of the paper and the results obtained. In the second section 
we start with describing the Poisson structure on T*GL(N)(z, z) that depends on two complex 
parameters k and a. Then we parametrize T*GL(N)(z, z) in a special way. The Poisson 
structure in new variables is ill defined due to the presence of singularities. To overcome this 
problem, we introduce an intermediate regularization. Removing the regularization we find 
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that only for a = 1/N the resulting Poisson structure is well defined. The value a = 1/N 
corresponds, in fact, to the case where only the SL(N)(z, z)-subgroup is centrally extended. 
The corresponding Poisson structure is described by two matrices r and r, which depend on 
N dynamical variables g,. It follows from the Jacoby identity (see Appendix A) that r is an 
iV-parameter elliptic solution to the classical Yang-Baxter equation (CYBE). It is worthwhile 
to mention that the main elliptic identities (see Appendix B) follow from the fulfillment of 
the CYBE for r. We expect that the matrix r is related to a special Frobenius subgroup in 



GL(N)(z,z) as it was in the finite-dimensional case [18[]. The Jacoby identity also implies a 
closed system of equations on r and f . 

We define a special matrix function L on T*GL(N)(z, z). We call this function the "L- 
operator" since as we show the Poisson algebra of L inherited from T*GL(N)(z, z) literally 



coincides with the one for the elliptic RS model || |20|. Thus, the classical L-operator algebra 



can be realized as a subalgebra of the algebra of functions on T*GL(N)(z,z). It turns out 
that the L-operator as a function on T*GL(N)(z, z) admits a factorization L = WP, where 
Pi = log Pi are the variables canonically conjugated to % and W belongs to some special 
subgroup in GL(N)(z, z). The Poisson bracket for the entries of W is given by the matrix r 
and coincides with the Sklyanin bracket defining the structure of the Poisson-Lie group. 

Although the quantum analogs of equations on r and f can be easily established, it is rather 
difficult to find the corresponding quantum R and R-matrices. The matter is that the matrices 
r and f have the complicated structure, r = r — s®I + I®s and r = f— s®/, due to the presence 
of the s-matrix. However, we observe that the classical L-operator algebra does not depend 
on s and, moreover, the matrices r and f also obey a closed system of equations. We show 
that this system arises as the compatibility condition for a new Poisson algebra. This algebra 
contains both the Poisson algebra of T*GL(N)(z, z) and the classical L-operator algebra as its 
subalgebras. In the third section, using this key observation, we pass to the quantization. 

We find the corresponding quantum R and .R-matrices as solutions to quantum analogs of 
the equations for r and f. In particular, the .R-matrix satisfies a novel triangle relation that 
differs from the standard quantum Yang-Baxter equation by shifting the spectral parameters 
in a special way. The Felder elliptic i? F -matrix naturally arises in our construction. It turns 
out that the i?-matrix is twist-equivalent to the i? F -matrix with the i?-matrix playing the role 
of the twist. 

Then we derive a new quadratic algebra satisfied by the "quantum" L-operator. This algebra 
is described by the quantum dynamical i?-matrices, namely, R, R F , and R: 

L42L1-R21L2 = L2-Ri2Li-R^. 

We show that the system of equations on R, R F , and .R-matrices arises as the compatibility 
condition for this algebra. We present the simplest representation of the quantum L-operator 
algebra corresponding to the elliptic RS model. We note that when performing a simple canon- 
ical transformation, the quantum L-operator coincides in essential with the classical L-operator 
found in [[K] . 



The quantum integrals of motion for the elliptic RS model were obtained in JT0[. In [gjj] , 
it was shown that any operator from the Ruijsenaars commuting family can be realized as the 
trace of a proper transfer matrix for the special L-operator that obeys the relation RLL = LLR 
with Belavin's elliptic i?-matrix |[26|| . We note that our L-operator is gauge-equivalent to L. 
It follows from this observation that the determinant formula for the commuting family pi] is 
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also valid for L. We show that any representation of our L-operator algebra is gauge equivalent 
to a representation of the relations RLL = LLR. 

In Conclusion we discuss some problems to be solved. 



2 Classical L-operator algebra 

2.1 Poisson structure of T*GL(N)(z, z) 

Let T T be a torus endowed with the standard complex structure and periods 1 and r. Denote 
by G a group of smooth mappings from T T into the group GL(N,C). Then g G G is a 
double-periodic matrix function g(z,z). The dual space to the Lie algebra of G is spanned by 
double-periodic functions A(z, z) with values in Mat (N, C). In what follows, we often use the 
concise notation g(z,z) = g(z) and A(z,z) = A(z). The group G admits central extensions 



G [ 22fl . The Poisson structure on T*G with fixed central charges reads 



{A 1 {z),A 2 {w)} = ~[C,A 1 (z)-A 2 (w)}6(z-w)-k(C-aI)-^5(z-w) (2.1) 

{ gi (z),g 2 (w)} = (2.2) 
{A 1 (z),g 2 (w)} = g 2 (w)C5(z-w), (2.3) 

where k, a are central charges and 5(z) is the two-dimensional ^-function. Here we use the 
standard tensor notation, and C is the permutation operator. 

One can consider the following Hamiltonian action of G on T*G 

A{z) -> T-\z)A(z)T(z) + kT' 1 {z)dT{z), (2.4) 
g(z) - T-\z)g(z)T(z). 

We restrict our consideration to the case of smooth elements A(z). Then generic element A(z) 
can be diagonalized by the transformation (|2.4j) [p3f : 



A(z) = T(z)DT- 1 (z) - kdT{z)T- 1 (z). (2.5) 

Here D is a constant diagonal matrix with entries D iy Di ^ Dj and T{z) is double-periodic. 
Matrix D is defined up to the action of the elliptic Weyl group. One can fix D by choosing the 
fundamental Weyl chamber. 

Matrix T(z) in eq. (|2.5| ) is not uniquely defined. Any element T(z,z) = T(z, z)h(z), where 
a diagonal matrix h(z) is an entire function of z, also satisfies ( J2.5| ). Demanding T(z, z) to be 
double-periodic, we obtain that h(z) is a constant matrix. We can remove this ambiguity by 
imposing the condition 

T(e)e = e, (2.6) 

where e is a vector such that e% = 1 Vi, and e is an arbitrary point on T T . In what follows, 
we denote the matrix T(z) that solves eq.( [2.5|) and satisfies eq.( [2.6|) by T e (z). Such matrices 
evidently form a group. 

Now we try to rewrite the Poisson structure (|2.1| ) in terms of variables T and D. Since 
Di are G-invariant functions, they belong to the center of ( |2.1| ) and, therefore, it is enough to 
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calculate the bracket {T e (z), T e (w)}. However, the straightforward calculation reveals that this 
bracket is ill defined. So, we begin with calculating the bracket {T £ (z), T v (w)}, where T e {z) 
and T v (w) satisfy ( |2.6|) at different points e and 77, 



m(z),T%(w)} = £ / rfVrfV l T } { ™] {A mn (z'),A ps (w')}. (2.7) 

mnps J 0/i mn {Z ) Ap S {W ) 



ST?, (z) 

To calculate the functional derivative SA 1,3 ^ , we consider the variation of (|2.5|) : 

= t(z)D - Dt(z) - fcdt(z) + d, (2.8) 

where X(z) = T^(z)5A(z)T(z), t(z) = T- 1 (z)5T(z) and d = 5D. 
First, from ( |2.8| ), we immediately obtain 

<)D ' 1 ^ <(z)T%z). (2.9) 



SA k i(z) (r - r 
Let us introduce the function $(z, s) of two complex variables 

s ) = ^^e^^e 2 ™^*. (2.10) 

Here a(z) and £(2) are the Weierstrass a- and ^-functions with periods equal to 1 and r. The 
function $(z, s) is the only double-periodic solution to the following equation: 

— 2lTtS 

d${z,s) + -&{z,s) = 2m5{z). (2.11) 

r — t 

It is also convenient to define $(z, 0) as follows: 

^z,0) = m($(z,e)-^=C(z)-2C(±)z + 2 ' Z ~ Z 
This function solves the equation 

2ni 

d<S>(z,0) = 2ni5{z) 



T — T 



T — T 



We introduce the notation = q-i — g 3 -, where % = |-4;Z)j. 

Using these functions, one can write the solution to ( |2.8| ) obeying the condition t(e)e = 

H : 

t(z) = d2 w(^(e - w, 7, ; ).V, ,(//•)/•;„ - $(> - w, q ij )X ij (w)E ij ). (2.12) 

Hereafter, we denote by «. 

Performing the variation of eq. fl2.12|) with respect to A mn {w) one gets 

5T e (z) ( 
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To compute the bracket (|2.7|) , one needs the following relation between T £ (z) and T v (z): 

T £ {z) =T ri (z)H r > £ , (2.13) 

where H V£ is a constant diagonal matrix. 
By direct computation, one finds 

Vf(z),T 2 >)} = Tt(z)T2{w){HFH?Wz,w) - affrw)), (2.14) 

K 



where 



r%(z,w) = X <1>! •//.'/,,)/•:„ /•-'/, • X <lM : ' ( 2 - 15 ) 
- 51 $(2 - 77, q^Eij ® E^j + ^ - e, ® £7 y (2.16) 



and 

/ £,? (z, w) = $(e — 77, 0) + $(cj - e, 0) + $(z - w, 0) - $(z — ry, 0). 

The bracket ( |2.14j ) has the r-matrix form with the r-matrix depending not only on coordinates 
qi but also on the additional variables H . 

In the limit rj — > e, one encounters the singularity. This shows that the variable T(z) is not 
a good candidate to describe the Poisson structure fl2.1|) . However, one can use the freedom to 
multiply T(z) by any functional of A. So, we introduce a new variable 

T £ (z) = T £ (z)(detT £ (e)f. (2.17) 

We use detT e (e) in the definition of T £ (z) in order to have the group structure for the new 
variables. 

Using the Poisson bracket (|2.14 ) one immediately finds 

~{Tf(z),T^)} = Tl(z)TUw)(HrHrr%(z,w)-ar(z,w) 

+ (31 ® tr 3 H £ ^HT4^e, w) + (3 tr s fff flfrgC*, 77) ® / 

+ 1 tr u Hl*H2 e rl2(e, 77) I ® /) , (2.18) 



since /(e, w) = /(z, 77) = 0. 

Now we are going to pass to the limit 77 — > eB For this purpose one should take into account 
the following behavior of H £V when 77 goes to e: H £V = 1 + (e — rj)h + o(e — 77), where h is 
a constant diagonal matrix being the functional of A. It turns out that there exists a unique 
choice for a and (3, namely, a = 1/N, (3 = —1/N, for which the singularities cancel and there 
is no contribution from the matrix h. In the limit 77 — > e = 0, i for these values of a and /3, one 
gets 

i{Tx(z), T 2 (w)} = T!(z)T 2 (u;)r 12 (z, «;). (2.19) 



1 Without loss of generality we assume that e and ij are real. 

2 The e-dependence can be easily restored by shifting both z and w by e. 
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Here the limiting r-matrix is given by 



where 



and 



r 12 (z, w) = r 12 (z, w) - s(z) <g) I + I <g> s(w) - ~ w, 0)1 <g> /, (2.20) 



r(z, w) = $(qij)Eu g> % + J] $(2; - iu, ® 

ij 

- £ «I>( :. ® /•-',; + 51 'I*!"'- ® E H (2-21) 

s W = ^ E (^(%)^i - $(^, <fo)3tf) • (2.22) 

ij 

Here we denote by the function the regular part of $(e, g^) for e — > 0: 

= 0. 

Note that both r and r are skew-symmetric: 1*12(2;, w) = —r2i(w,z). 

A natural conjecture is that the r-matrix obtained satisfies the classical Yang-Baxter equa- 
tion 

[[r, r]] = [r 12 (z h z 2 ), r X3 {z l: z 3 ) + r 23 (z 2 , z 3 )] + [r^fa, 2:3)^23(^2, z 3 )] = 0. (2.23) 

It can be verified either by direct calculation or by considering the limiting case of the Jacoby 
identity for the bracket (|2.18|) as is done in the Appendix A. Thereby, the r-matrix ( 2.20|) is an 
iV-parameter solution of the classical Yang-Baxter equation. 

Let us note that, as one could expect, the condition det T(0) = 1 is compatible with the 
bracket ( |2.19|) , since detT(;z) is a central element of algebra ( [2.19|) . 

Remark that the choice a = 1/N corresponds to the case where only the sl(N)(z, z)- 
subalgebra is centrally extended. In terms of T(z) ((3 = —1/N), the boundary condition looks 
like T(0)e = Ae and det T(0) = 1. One can also check that the field A(z) defined by ( [2.5|) with 
the substitution T(z) for T(z) obeys Poisson algebra (|2.1| ). 

The next step is to consider the special parameterization for the field g(z). To this end, we 
introduce A(z): 

A = gAg- 1 - kdgg' 1 + A tr dgg' 1 . (2.24) 
One can check that A(z) Poisson commutes with A(w) and obeys the Poisson algebra: 

{A 1 (z),A 2 (w)} = -~[C,A 1 (z)-A 2 (w)]6(z-w) + k(C-±I)-^S(z-w) 

{Ai(z),^H} = Cg 2 (w)5(z-w). (2.25) 
Now we factorize A(z) in the same manner as it was done for A(z), 

A(z) = \]{z)D\J-\z) - kd\]{z)\]-\z), (2.26) 
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where U(z) satisfies the boundary condition U(0)e = Ae and detU(O) = 1. Obviously, V(z) 
Poisson commutes with T(w) and satisfies the Poisson algebra 

-{Ui(«), V 2 (w)} = -V 1 (z)V 2 (w)r 12 (z, w). (2.27) 

K 

One can find from fl2.24| ) and (|2.26 ) the representation for the field g, 

g(z) = (detsO^U^PT- 1 ^), (2.28) 

where P is a constant diagonal matrix. 

Computing the determinants of the both sides of eq. (|2.28|) one gets 

detP = det (T(z)/U(*)). 

Since the l.h.s. does not depend on z and det T(0) = detU(O) = 1 we obtain that det P = 1 
and detT(z) = detU(z). 

Calculating the Poisson brackets of P with P and Q = diag(gi, . . . q^) in the same manner 
as above one reveals that 

{Pi,P 2 } = 0, (2.29) 
-{Qi,P 2 } = P 2 {J2E ii ®E ii -^-I®I). (2.30) 

n 

In fact, it means that logP; = Pi — jjJ2iPii where Pi are canonically conjugated to q it 
For the remaining Poisson brackets of P with the fields T, U, we have 

-{Ti(*),P 2 } = Ti^Parxa^), (2.31) 

K 

-{U!(z),P 2 } = U!^)P 2 r 12 (z). (2.32) 

K 

Here 

fi 2 {z) = f 12 {z) - s{z) <g> I - -I ® Yl 'H'/,,!/^,- (2-33) 

ij 

where we introduced the f-matrix: 

fi 2 (z) = ]T §{q l3 )E u ® E J3 - ]T $(z, ® %• (2.34) 

ij ij 

To complete the description of the classical Poisson structure of the cotangent bundle we 
present the Poisson bracket of det g with other variables: 

-{Q, det g(w)} = detg(w), {P, det g(w)} = 0, 

K 

-{T( z), det g(w)} = - det g(w)T(z)(<$>(z - w, 0) + 0)), 

K 

~{U (z), det g(w)} = - debg(w)U(z)($(z - w, 0) + 0)). 



Recall that the Jacoby identity for bracket (|2.19 ) reduces to the classical Yang-Baxter 



equation for the r-matrix. As to the Poisson relations ( |2.31| ) and fl2.32| ), one finds that the 



Jacoby identity is equivalent to the following quadratic in r equation: 

[f l2 {z), f 13 (z)} - P?{f n {z), P 3 } + P 2 _1 {f is (z), P 2 } = (2.35) 
and the equation involving r and f , 

[r 12 {z,w),f 13 (z) + f 23 {w)) + [f 13 {z),f 23 {w))-P 3 l {f 12 {z,w),P 3 } = 0. (2.36) 



One can check by direct calculations that the matrices r and r given by Q2.2UQ and fl2.33| ) do 
solve these equations. 



Finally, we remark that the fields A(z) and g(z) defined by ( |2.5|) and (|2.28| ) obey Poisson 
relations 



Let us note that the Poisson relation for the generator W(z) = T 1 (z)\J(z) turns out to 
be the Sklyanin bracket: 

-{Wi(*),W 2 (tu)} = [r 12 (z,w),W 1 (z)W 2 (w)\, (2.37) 

which, therefore, defines the structure of a Poisson-Lie group. This group is an infinite- 
dimensional analog of the Frobenius group appeared in [18| , where the Poisson-Lie group struc- 
ture was related to the existence of a non-degenerate two-cocycle on the corresponding Lie 
algebra. It would be interesting to find a similar interpretation in the infinite-dimensional case. 

2.2 Classical L-operator 

In this subsection, we define a special function on the cotangent bundle, which we call the 
classical L-operator. The motivation to treat this function as the L-operator is that the Poisson 
algebra of L is equivalent to the one found in || for the L-operator of the elliptic RS model. 
Denote by L the following function 

L(z) = T-\z)g(z)T(z) = (det g(z))^T-\z)V(z)P. (2.38) 

By using the formulas of the previous subsection, one can easily derive the Poisson bracket of 
L and Q: 

-{Q u L 2 (z)} = L 2 (z) ]T E u (2) E ti , (2.39) 
K i 

and the Poisson algebra of the L-operator, 

-{L 1 (z),L 2 (w)} = r 12 (^,«;)L 1 (z)L 2 («;) (2.40) 

K 

+ L 1 (z)L 2 (w)(f 12 (z) - f 2X (w) - r X2 {z,w)) 
+ Li(z)r 2 i(w)L 2 (w) - L 2 (w)f 12 (z)Li(z). 

Clearly, the generators Q and L form a Poisson subalgebra in the Poisson algebra of the 
cotangent bundle. An important feature of this subalgebra is that I n {z) = trL n (z) form a set 
of mutually commuting variables. 
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Just as in the finite-dimensional case ]18| , one can see from (|2.39| ) that the L-operator admits 
the following factorization: L(z) = W(z)P, where Q and logP are canonically conjugated 
variables, the W-algebra coincides with ( |2.37|) , and the bracket of W and P is 

-{Wx(z), P 2 } = -P 2 [ria(«), W x {z)\. (2.41) 

K 

In fact, everything what we need to quantize the L-operator algebra ([2.40 ) is prepared. The 
problem of quantization is reduced to finding the quantum R and R-matrices satisfying the 
quantum analogs of eqs.( |2.23D , ( |2.35| ), and ( |2.36| ), 

R 12 Czi,* 2 )R2i(*2,*i) = 1, (2.42) 

Rn{Zl, Z^jH^Zl, Z 3 )K 23 (z 2 , Z 3 ) = R-23(22,^3)Rl3(2l, 23)Rl2(2l,22), (2.43) 

R 12 (^,^)R 13 (zi)R 23 (2 2 ) = R 23 (z 2 )Ri3(2i)P 3 - 1 Ri 2 (zi,^)P3, (2.44) 
R 12 (^)P 2 - 1 R 13 (^)P 2 = R 13 (z)P 3 - 1 R 12 (^)P 3 . (2.45) 

These matrices are assumed to have the standard behavior near H = 0: 

R= l + hr + o(h), R=l + hf + o(h), 

where h is a quantization parameter. 

The problem formulated seems to be rather complicated due to the presence of the s-matrix 
in the classical r and f-matrices. However, Poisson algebra ( |2.4Q|) possesses an important 
property allowing one to avoid the problem at hand. Namely, the matrix s(z) coming both in 
r and r drops out from the r.h.s. of ( |2.4(J| ). Thereby, eq. fl2.40p can be eventually rewritten as 

~{Li(z), L 2 (w)} = r^z.w^^L^w) (2.46) 

- L 1 (z)L 2 (w)(ri 2 (z,w) + T2i{w) - n 2 (z)) 
+ L l (z)f 2l (w)L 2 (w) - L 2 (w)f 12 (z)L l (z). 

Moreover, if we denote by rf 2 the sum 

rf 2 (z,w) = r l2 (z,w) + f 21 (w) - f 12 (z) , (2.47) 

then using (|2.21|) and (|2.34| ) we obtain 

r{ 2 (z -w) = -J2 ®(Qij)Eii ® Ejj + $(« - w, q^E^ <g> E jt . (2.48) 

ij ij 

In this expression one can recognize the elliptic solution to the classical Gervais-Neveu-Felder 
equation |rj| |15[ : 

[rf 2 (zi - z 2 ),r( 3 (z! - z 3 ) + r^{z 2 - z 3 )} + [rf 3 (*i - z 3 ), rf 3 (z 2 - z 3 )} (2.49) 
-P 3 - 1 {rf 2 (^ 1 - z 2 ),P 3 } + P 2 - 1 K 3 (2i - z 3 ),P 2 } - Pr 1 {^ 3 (22 - 2 3 ),Pi} = 0. 

In fact, r F emerges as the semiclassical limit of the quantum R- matrix found in fli"5f| . 

The absence of the s-matrix in the resulting L-operator algebra and the appearance of 
the r F -matrix show that there may exist a closed system of equations involving only r- and 



10 



f-matrices in the classical case, and R- and .R-matrices in the quantum one. In the next 
subsection we find the desired system of equations and describe a Poisson structure for which 
these equations ensure the fulfillment of the Jacoby identity. 

Note that the algebra ( |2.46| ) literally coincides with the one obtained in [2"U| by using of the 



Hamiltonian reduction procedure. A mere similarity transformation of L turns algebra ( |2.46| ) to 
the one previously found in ||. In contrast to || where fl2.46p was derived by direct calculation 



with the usage of the particular form of the L-operator for the RS model, our treatment does 
not appeal to the particular form of L. 

2.3 Quadratic Poisson algebra with derivatives 



In the first subsection, we obtained the matrices r and r obeying system of equations (2.23) 



( |2.35| ) and (|2.36|) . Clearly, these equations are not satisfied when substituting r and r for r 



and r. However, computing the l.h.s. of these equations after this substitution we arrive at 
surprisingly simple result: 

[rnizi, z 2 ), ri 3 (2i, z 3 ) + r 23 (z 2 , z 3 )) + [rjz(zu z 3 ), r 23 (z 2 , z 3 )) = (2.50) 
-(<9i + d 2 )r l2 (z 1 , z 2 ) + (<9i + d 3 )n 3 (zi, z 3 ) - (d 2 + d 3 )r 23 (z 2 , z 3 ), 

[fMz), f 13 (z)} - P^in^z), P 3 } + P^in^z), P 2 } = -d{f 12 (z) - f 13 {z)), (2.51) 

and 

[rv»(zi, z 2 ), f 13 (zi) + f 23 (z 2 )] + [fi 3 (zi), f 23 (z 2 )] - P^in^Z!, z 2 ), P 3 } 

= -{d x + d 2 )r 12 (z l7 z 2 ) + d x f l3 (zi) - d 2 f 23 (z 2 ). (2.52) 

Here d — where x = Re z. Note that eqs.( ^.35|) and ( ^.36| ) are formulated with the help 
of P. However, since all the matrices depend only on the difference qij = qi — qj, we simply 
replace P by P. 

Comparing eqs.( p.50"| - p.52| ) for r and f with (|2.23| ), ( |2.35j ), and ( p.36|) for r and f , we come to 



the conclusion that the s(,z)-Hiatrix coming in r and f effectively plays the role of the derivative 
with respect to the spectral parameter. 

It is worth mentioning that eqs.( p.50"|) -( p.52|) obtained for r and f can be rewritten in the 



same form as eqs. (|2.23|) , (|2.35|) , and (|2.36| ) if we replace r and r by r i2 — d± + d 2 and r\ 2 — d\. 



In particular, for (|2.50| ), we have 

[ri2 ~di + d 2 , r 13 -di+ d 3 ] + [r 13 -dt + d 3 , r 23 -d 2 + d 3 ] + [r 12 -di+d 2 , r 23 -d 2 + d 3 ] = 0. (2.53) 

Thus, ri2 — d\ + d 2 is a matrix first-order differential operator satisfying the standard classical 
Yang-Baxter equation. Using this fact we write down the Poisson algebra generated by the 
fields T(z), U(z), Q and P, having eqs.( p.50|) -( p.52| ) as the consistency conditions: 

-flK^TsW} = T 1 (z)T 2 (w)r 12 (z,w)+T[T 2 -T 1 T^ (2.54) 

K 

-{Ux(z),U 2 (w)} = -U 1 {z)U 2 {w)r l2 {z,w)-U[U 2 + U 1 U' 2 , (2.55) 
-{T x {z),P 2 } = P 2 T 1 (z)f 12 (^) + P 2 T 1 / (^), (2.56) 

K 
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;{U 1 (Z),P 2 } = P 2 U 1 (z)f 12 (z)+P 2 U[(z), 

{Q U P 2 } = P 2 Y,E U ®E U , {Q 1 ,T 2 } = {Q 1 ,U 2 } = 0, 



(2.57) 
(2.58) 



where T' = dT. 

It is worth mentioning that the Poisson structure (|2 . 54j - p~58|) is not compatible with the 
boundary condition T(0)e = Ae. 

Let us note that there exists a Poisson subalgebra of Poisson algebra ( |2.54| - f2~5"8| ), formed by 
the generators: 

A{z) = T{z)DT-\z) - kdT{z)T'\z), g(z) = U(z)PT~ 1 {z) 

that coincides with the Poisson algebra of the cotangent bundle with the central charge a = 0. 

Defining the L-operator as L(z) = T~ 1 (z)g(z)T(z) = T~ 1 (z)U(z)P, we get for L algebra 
( |2.46| ) obtained previously. As in the previous subsection the commutativity of I n (z) follows 
again from the one of g(z). 



The main advantage of Poisson algebra (|2.54-|2.58|) is that it can be easily quantized. 



3 Quantization 



3.1 Quantum L?-mat rices 

we 



In this section, following the ideology of the Quantum Inverse Scattering Method p4| , |25 
quantize the classical r and f-matrices and derive the quantum L-operator algebra. 

We start with quantization of the relations ( 2.50|) -( 2.52 ). Let T(z), U(z) be matrix gener- 
ating functions being the formal Fourier series in variables x and y: 

T(z) = J2 T mn e 2 ^ mx+n ^ , U(z) = J2 U mn e 2m( - mx+n ^ , 

mn mn 

where z = x + ry. 

Denote by A a free associative unital algebra over the field C generated by matrix elements 

of the Fourier modes of T(z), U (z), and by the entries of the diagonal matrices P and Q modulo 
the relations 

T 1 (z)T 2 (w-h)=T 2 (w)T 1 (z-h)Rx2(-h,z,w), (3.1) 

U^U^w + h) = U 2 (w)U!(z + h)R 12 {h, z, w), (3.2) 

T x {z + h)P 2 R 12 (h, z) = P 2 Ti(z), (3.3) 

U x (z + h)P 2 R 12 (h, z) = P 2 U x {z), (3.4) 

[Q 1 ,P 2 \ = -hP 2 '£E ii ®E ii , (3.5) 

i 

[T l (z),U 2 (w)] = [T^Qz] = [U x (z),Q 2 ] = [P U P 2 ] = [Q U Q 2 ] = 0. 

Here R(h,z,w) and Ri 2 (H,z) are double-periodic matrix functions of spectral parameters. 
These functions also depend on the coordinates qi and have the following semiclassical be- 
havior at h = 0: 

R=l + hr + o(h), R=l + hf + o(h). (3.6) 
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The next step is to find the conditions on R and R that ensure the consistency of the defining 
relations for A. In the sequel we often use R(z, w) as a shorthand notation for R(h, z, w). 

First, we write down the compatibility condition for algebra ( |3.1| ) or ( |3.2| ), which reduces 
to the Quantum Yang-Baxter equation with spectral parameters shifted by h, 

R X2 (z, w)R 13 (z -h,s- h)R 23 (w, s) = R 2 3(w — h,s — h)R 13 (z, s)R X2 {z — h,w — %). (3.7) 

Analogously to the classical case, one can introduce the following matrix differential operator 
lZ(z,w) = e h 9^ R(z,w)e~ h 8^ in terms of which eq.( |3.7|) reads as the standard Quantum Yang- 
Baxter equation 

K 12 (z, w)TZ 13 (z, s)TZ 23 (w, s) = 1Z 2 3(w, s)TZx 3 (z, s)K 12 (z, w). (3.8) 
Relation ( |3.1| ) also requires the fulfillment of the "unitarity" condition for R, 

Ri 2 (z,w)R 21 (w,z) = 1. (3.9) 
Analogously, we find the following compatibility conditions for (|3.3j ): 

P.i l R 12 (z)P 3 Ry3(z -h)= P 2 - l R 13 (z)P 2 R 12 (z - K) (3.10) 

and 

P^Ruiz, w)P 3 R 13 (z - h)R 23 (w) = R 23 (w - h)R 13 {z)R 12 {z -h,w-h). (3.11) 
Now taking into account ( |3.6| ) one can easily see that in the semiclassical limit 

--{•,•} = Jim 

relations ( |3TT| - |3.5| ) determine Poisson structure ( p.54| - |2.5^ ), while eqs.( p.7| ), ( |3.10| ), and ( |3.11| ) 
turn into ( p.50| ), ( p.51|) , and (p.52|) , respectively, in order h 2 . In the first order in h, the unitarity 
condition ( |3.9| ) requires r to be skew-symmetric. Hence, the algebra A with defining relations 
( j3TTD -( j3T5D , where R and R are the solutions of ( j3T7| - [3.11| ) obeying (p.6|) , is a quantization of the 
Poisson structure ( |2.54| - |2.58"D . 

Now we are in a position to find the matrices R and R explicitly. We start with the i?-matrix 
for which we assume the following natural ansatz: 

fR(h, z,w) = J2 $(^1, Qij + ^2)En ® Ejj + J2 H z ~ w + ^3, Qij + fri)Eij ® Eji 

ij ij 

- ®(z + h 5 , qij + h 6 )E i:i <g> Ejj + ^ $(w + h 7 , g^- + h 8 )E jj ® E {j . (3.12) 

ij ij 

This form is compatible with the structure of the classical r- matrix. Here Hi, . . . ,hg are ar- 
bitrary parameters that should be specified by eqs.( |3.7| ) and ( |3.9| ), and / is a scalar function 
that may depend only on hi and spectral parameters. It turns out that the parameters hi are 
almost uniquely fixed by the unitarity condition (|3.9|). Substituting ( |3.12| ) into ( |3.9| ) and using 
the elliptic function identities we obtain 

h 2 = h 3 = hi = h 6 = hs = o, hs = hx + h 7 , f 2 (z, w) = V(hi) - v(z - w), 
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where V(z) is the Weierstrass P-function. Now it is a matter of direct calculation to check that 
eq.( |3.7| ) holds for hi = h. The remaining parameter h-j is inessential since it corresponds to an 
arbitrary common shift of the spectral parameters z and w. In the sequel, we choose T17 = 0. 
Therefore, the obtained solution to fl3.7| ) and ( |3.(j| ) reads as follows: 



f(z,w)R(H,z,w) = ^$(h,q ij )E ii ®E jj + Y / $(z-w,q ij )E ij (g)E^ 

ij ij 
- ®( Z + h ' 1ij) E ij ® E jj + 1ij) E jj ® E ij: 



(3.13) 



'J 



where f(z, w) = JV(ti) — V{z — w). One must be careful in the definition of R(—h, z, w). This 
matrix is defined by ( |3.13|) with the replacement H — > — 7i and / — > — /. Therefore, R(H,z,w) 
and R(—H,z,w) are related as 



Ri2{—h, z, w) = R2l(fo, w — h, Z — h). 
To find the i?-matrix, we adopt the following ansatz: 



(3.14) 



1 ^ 



#12 (ft, z) = 2 $ (^> ®i + fta)-B« ® - J] ®( z + ha, q-ij + ^4 + ^5)-^ <g> (3.15) 



It has almost the same matrix structure as the classical r- matrix. Since eq.( |3.10| ) is easier to 
deal with than eq.( |3.11|) , we first substitute (|3.15|) into eq. (|3.10|) thus obtaining R: 

-^rR 12 (h, z) = ®(h, Qij)Eii ® e^ - + h 3 , <i, J )E, J ® % - $(z + ft 3 , -K) E u ® E u- 

(3.16) 

where /13 remains unfixed. 

Eq.( |3.11| ) involves both the R- and .R-matrices and is independent on (|3.7| ) and (|3.10| ). One 
can verify by direct calculations that R and R given by eqs. fl3.13| ) and ( |3 . 1 6| ) also satisfy ( |3.11| ) 
as soon as = h. 

One can easily check that in the case of real h, the matrices R and R have the proper 
semiclassical behavior (|3.6|) . 

In what follows we also need the R -matrix, 



-TTxR^ih, z) = -J2 $ (-ft> Qij + n ) E n ® E jj + E $ (z, Qij + h)Eij ® E. 

(T \' 1 ) H ii 



(3.17) 



It would be of interest to mention that just as in the rational case without the spectral 
parameter [IS|, one can introduce the formal variable W(z) = T^ 1 (z)U(z) with permutation 
relations following from ( |3.1| - |3.5| ): 



R 12 (z,w)W 1 (z)W 2 (w + h) 
W 1 (z + h)P 2 R 12 (z) 



W 2 (w)W 1 (z + h)R 12 (z,w) 

P2Rl2(z)W 1 ( Z ). 



(3.18) 
(3.19) 



In analogy with the rational case, it is natural to treat eq.( |3.18 ) as the defining relation of the 
quantum elliptic Frobenius group. 
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3.2 Quantum L-operator algebra 

Just as in the classical case, we introduce a new variable: 

L(z) = T-\z)U(z)P = W(z)P, (3.20) 

which we call a quantum L-operator. Using the relations of the algebra A one can formally 
derive the following algebraic relations satisfied by the quantum L-operator: 

[Qi,L 2 (z)} = -%L 2 {z)Y,Eii®E iU (3.21) 

i 

Ri 2 (z,w)Li(z)R 2 i(w)L 2 (w) = L 2 {w)Ri2{z)Li(z)R 2 i{w - h)R 12 (z -h,w — h)R 12 (z - h). 

(3.22) 



In spite of the fact that L has the form L(z) = W(z)P, we can not reconstruct from eqs. (|3.21 



and ( |3.22| ) the relations ( |3.18| ) and ( |3.19| ) for W and P. So, in the sequel, we do not assume 
any relations on W and P. 
Let us define 

Rf 2 {z,w) = R 2l {w)R 12 {z,w)R l2 \z). (3.23) 

Then, by using the explicit form of the R- and i?-matrices and elliptic function identities, we 
obtain 

fRf 2 (z -w) = -J2 qij)Eu <g> Ejj + w, Qi^Eij <g> Eji + <£>(z -w,h)J2 E u ® E„, 

(3.24) 

which is nothing but the i?-matrix by Felder [I5]j, i.e., an elliptic solution to the quantum 
Gervais-Neveu-Felder equation []16l p~5f] : 



P^R^w - s)P 1 Rf 3 {z - s)P.^RUz - w)P 3 = Rf 2 {z - w)P 2 l RU* ~ s)P 2 i^(^ - s). (3.25) 
Recall that one feature of R F is the "weight zero" condition: 

[P 1 P 2 ,Rf 2 (z-w)}=0. (3.26) 



Developing R F in powers of h, we have R F = 1 + hr F + o(h), where r F is given by ( [2.48|) . 

Let us stress that in our consideration R F arises to account for the explicit form of R and R, 
and that the Gervais-Neveu-Felder equation does not follow from system ( |3.7| - [3.11 ). Formula 
( [3.231 ) shows that the matrix R plays the role of the twist, which transforms the matrix R(z, w) 
- a particular solution of (|3.7| ) - into such solution of ( |3.25| ) which depends only on the difference 
z — w. 

Thus, the quantum L-operator algebra ( p.22[ ) can be presented in the following form: 

J R 12 (2,w;)L 1 (^)L 5 21 («;)L 2 (w;) = L 2 (w)R 12 (z)L 1 (z)R F 2 (z - w). (3.27) 

The quantum L-operator algebra seems to be automatically compatible as A is compatible. 
However, a simple analysis shows that A and the algebra (|3.27|) admit different supplies of 
representations. In particular, the simplest representation for L we present below does not 
realize the algebra ( |3.18|) , ( |3.19| ). Therefore, we find it necessary to give a direct proof of the 
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compatibility of fl3.27Q . In this way, we come across eq.( |3.25|) and discover a new relation 
involving R F and R. To this end, let us multiply both sides of ( |3.27| ) by 



P 2 l R 3l {s + h)P 2 R 32 {s)L 3 {s) 

and subsequently using eq. (|3.27|) we transform the string Li ■ ■ 
For the l.h.s., we have 

J R 12 (^,w)L 1 ^ 21 (w;)L 2 P 2 - 1 E3 1 ( S + h)P 2 R 32 {s)L 3 = 
R 12 (z, w)L 1 R 21 (w)R 31 (s + K)L 2 R 32 (s)L 3 = 
R 12 (z,w)L 1 R 21 (w)R 31 (s + h)R 32 (s,w)L 3 R 23 (w)L 2 R 23 



■ L 2 ■ ■ ■ L 3 into L 3 ■ ■ ■ L 2 ■ ■ ■ L\. 



w 



R 12 (z,w)R 32 (s + k,w + ^)L 1 ^3i(s)L 3 P3 _1 ^2i(w + h)P 3 R 23 (w)L 2 R, 
R 12 (z,w)R 32 (s + h,w + h)R 31 (s,z)L 3 R 13 (z)L 1 x 
Rf 3 (z - s)P. 3 ~ 1 R 21 (w + h)P 3 R 23 (w)L 2 R F 3 (w - s). 



23 



(3.28) 



At this point, we interrupt the chain of calculations by remarking that the next step implies 
the possibility to push R F somehow through P^R^iw + h)P 3 R 23 (w). It can be done by virtue 
of the following new relation involving R F and R: 



Ru(z)P 2 1 R 31 (w)P 2 R 32 (w -K) = Pr 1 Rz 2 (w)P 1 R 31 (w - h)R 



121 



(3.29) 



which can be checked directly by using the explicit forms (|3.24|) and (|3.16|) of R F and R 
respectively. 

Now we pursue calculation (|3.28|) with the relation ( p.29| ) at hand. 

^12(^,^)^32(3 + h,w + h)R 31 (s, z)L 3 R 13 {z)R 23 {w + h) x 
UR^w^P^R^z - s)P 2 R F 3 {w - s). 

As to the r.h.s., the same technique yields 



L 2 R l2 {z)LxR 
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w)P 2 - 1 R 31 {s + h)P 2 R 32 (s)L 3 



L 2 R 12 {z)R 32 {s + %)L{R 3 i(s)L 3 P 3 1 R F 2 (z - w)P 3 = 

L 2 R 12 (z)R 32 (s + h)R 31 (s, z)L 3 R 13 {z)L x R F 3 {z - s)P 3 l R F 2 (z - w)P 3 

= R 31 (s + h,z + h^R^L^R^z + h)P 3 x 
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s)P 3 l R F l2 {z-w)P 3 



R 3 i(s + h,z + h)R 32 (s, w)L 3 R 23 (w)R 13 (z + h)L 2 R 1 2(z)L 1 x 

P^R^w - s)P 1 R F 3 (z - s)P 3 l R F l2 (z -w)P 3 = 
#31 (s + h,z + h)R 32 (s, w)R 12 (z + h,w + h)L 3 R 13 (z)R 23 (w + K)Li x 



(3.30) 



Rniw^Ri^w - z)P^ l R F 3 {w - s)P 1 R F 3 {z - s)P 3 l R F 2 {z - w)P 3 . 
Therefore, comparing the resulting expressions we conclude that the compatibility condition 



for the L-operator algebra ( [3.271 ) reduces to four equations (|3.7| ), (|3.11| ), Q3.24Q , and (|3.29| ). 

The existence of the Poisson commuting functions I n {z) in the classical case implies that the 
commuting family should exist in the quantum case as well. It should be the intrinsic property 
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of the algebra ( 3.27 ) itself, without referring to the explicit form of its representations. Let us 
demonstrate the commutativity of the simplest quantities txL{z) and trL -1 (;z) postponing 
the discussion of the general case to the next section. To this end, we need one more relation 
involving the matrices R F , R and R. 

In analogy with the rational case, it is useful to introduce the variable g(z) = U(z)PT~ 1 (z). 
Calculation of the commutator [gi(z) , g 2 (w)] with the help of the defining relations of A results 
in 



U 2 (w)U 1 (z + h) (Ri 2 {h, z, w)P l R 21 {w)P 2 R l2 {z - h)Ri 2 (-h, z, w) 
P 2 R l2 {z)P{R 21 {w - %j) T 2 \w - K)T{\z). 



[gi(z),g 2 (w)) 

When the spectral parameter is absent, the algebra A allows one to establish a connection 
with the quantum cotangent bundle (see for details). Then, in particular, the quantity 
[gi,g 2 ] is equal to zero. In the case at hand, we can not construct a subalgebra of A that is 
isomorphic to the quantum cotangent bundle. However, one can note that in the elliptic case, 
the commutativity of g(z) with g(w) follows from the identity 

R l2 {h, z, w)P 1 R 21 {w)P 2 R 12 {z - h)R 12 (-h, z, w) = P 2 R 12 (z)P 1 R 21 (w - K). 

Using the definition of R F , eq.( [3.14 ), and the "weight zero" condition ( 3.26 ), the last formula 
can be written in the following elegant form 



R 12 (z,w) = P 2 R l2 {z)P 2 l R( 2 {z - w)P 1 R 21 HPr 1 . 



(3.31) 



Identity (|3.31|) plays the primary role in proving the commutativity of the family trL(^) . To 
prove the commutativity, let us multiply both sides of 

L 2 (w)Ri 2 (z)L 1 (z) = Ri 2 (z,w)L 1 (z)R 21 (w)L 2 (w)R F Jw - z). 



by P 2 R 12 (z)P 2 and take the trace in the first and the second matrix spaces. We get 



tr 12 P 2 P 12 1 (^)P 2 ^ 1 L 2 (w;)P 12 (^L 1 (^ = 
tr u P 2 R~ 2 \z)P 2 1 R 12 (z,w)L 1 (z)R 21 (w)L 2 (w)R F 1 (w - z) 



It is useful to write P 12 in the factorized form 



(3.32) 



(3.33) 



where 



-a(h)<$>(-%, qij + %)E U + cr(h)<£>(z, + h)E i:j . 
Then the l.h.s. of (13.321) reads as 



r n {.PjRij Pj 1 



v 
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Taking into account that R\ 2 is diagonal in the second matrix space and using the cyclic 
property of the trace, we obtain 



tr ,2! /';/.',/ ; ® /)(/ ® LEjj)Ri 2 (z)Li 



Since L = VVP, where all entries of W commutes with g i; we arrive at 

J2 tr : AI , ~I!;/1' J ' ® /)(/ ® WP,P^ #12(^1 



£tr 12 L 2 (p: 



%j)P 



L2 



z)Li = trL(w) trL(z) 



As to the r.h.s. of (|3.32| ), we use identity (|3.31|) to rewrite it as 



tr i 2 R$ 2 (z - w)P l R 2 l{w)P{ 1 L l {z)R 2l {w)L 2 {w)R 



21 



w 



Having in mind that R 2 i is diagonal in the first matrix space and taking into account the 
property ( |3.26| ) one can easily see that under the trace sign, the matrix R F can be pushed to 
the right where it cancels with R 21 . Therefore, we get 



ti 12 P 1 R 21 \w)P 1 - 1 L 1 (z)R 21 (w)L 2 (w). 



(3.34) 



Now applying to this expression the technique we used above for the l.h.s of ( |3.32| ) we con- 
clude that eq. (|3.34|) is equal to trL(z) tr L(u>) . Thus, we proved that tr L(z) commutes with 
tr L(w) . Quite analogously one can prove that trL _1 (^) commutes with trL _1 (w) and with 
tr L{w) . 

Now we give an example of the simplest representation of algebra (|3.21 ) and (|3.27 ) associated 
with the elliptic RS model. Namely, the following L-operator satisfies algebra ( |3.27|) : 



L(z) = J2 ®( z > % + l)h P i E iv 



(3.35) 



where 

&; = II*(7 5 <to)- (3-36) 

Here the parameter 7 is a coupling constant of the elliptic RS model. This can be checked 
by straightforward calculations. Some comments are in order. The L-operator of the form 
( |3.35| ) was already appeared in []20| as a result of the Hamiltonian reduction procedure applied 
to T*GL(N)(z, z). To guess the explicit form of bj one should note that in the rational limit 
algebra (|3.27| ) tends to the one obtained in fTS] , where the coefficients bj are found to be 



n 



Qaj + 7 



(3.37) 



Therefore, it is natural to assume that the elliptic analog of (|3.37|) is given by ( |3.36 ). 

It is worthwhile to mention that bj are not uniquely defined since one can perform a canonical 
transformation of (q, p)-variables. In particular, the variables 



r = -pr g(gaj + 7) 

3 Ih^i)^) 



(3.38) 
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are related to bj by the canonical transformation — > and 

where a = 2C(l/2) 7 - z£j. 

We call this L the quantum L-operator of the elliptic RS model. Indeed, taking the Hamil- 
tonian to be H = trL(z) one can see that the quantum canonical transformation of the form: 

where P t R is the momentum in the Ruijsenaars Hamiltonian, turn H into the first integral 
Si from the Ruijsenaars commuting family [JLO]. Moreover, after the canonical transformation 
( p.39| ), the L-operator Q3.35| ) coincides in essential with the classical L-operator of the RS model. 
The generating function for the commuting family in terms of L can be written as 

/(*, /i) =: det {L{z) - n) := £(-/*)"-***(*), ( 3 - 40 ) 

k 

where the normal ordering :: means that all momentum operators are pushed to the right. 
It follows from the results obtained in the next section. 



4 Connection to the fundamental relation RLL = LLR 

In this section we establish a connection of the quantum L-operator algebra ( |3.27| ) with the 
fundamental relation RLL = LLR. 

In the operators from the Ruijsenaars commuting family were obtained by using a 
special representation L of the algebra 



RfJz - w)L 1 (z)L 2 (w) = L 2 (w)L 1 (z)R? 2 (z - w 



(4.1) 



where R B (z) is Belavin's .R-matrix being an elliptic solution to the quantum Yang-Baxter 
equation ||26|| . The explicit form of R B we use here can be found in |27|]. For reader's convenience 
we recall a construction of L [^R 

Denote by h* the weight space for sZjv(C) that can be realized in C N with a basis e«, 
< €i,€j >= Sij, as the orthogonal complement to Y,iLi e i- Let be the orthogonal projection 

£ 1 V^JV 

01 Cfc! 6k — tk — jj 2^i=l e i- 



For each q e h* one can introduce the intertwining vectors fl32| , |33 



where 



^■(tt- < Q^k >)/iv(r) 



' N 



(4.2) 



exp 27rz 



ne#-j+jvz 



. 1 n 2 

n z H — H r 

V 2 ; 2iV 



and 77 ( 



r 



P 



1/24 IT =ll 



1 — p m ) is the Dedekind eta function with p = exp 2mT. 
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(4.3) 



Following |28| we denote by <fi(z) q q +h!Lk i the entries of the matrix inverse to 0(z)« +nE Y Then the 
orthogonality relations read as follows 

n 

In the sequel, the following formula [2~I] will be of intensive use 



fc=i 



m=l 



z+ < q', €j > - < q, e { >) -j-j- 0(< g', e r > - < q, €j >) ^ ^ 



0(< g', >-<q', €j > 
Here #(z) denotes the Jacoby ^-function 

_ _^ e 27ri( 2 +|)(n+i)+i7rr(n+i) 2 = 0'(O)e~ C( ^ 2 (7 



3 . 



It is shown in 29, 30] that the L-operator 



v 



(4.5) 



k=l 



acting on the space of functions on h* satisfies relation fl4.1| ). This L is an N x N generalization 



of the 2x2 Sklyanin L-operator [31 



The intertwining vectors <fi(z) q q +hek j coming in the definition of L relate the matrix R B with 

the Boltzmann weights for the face model. Recall that the nonzero Boltzmann weights 
depending on the spectral parameter z are explicitly given by 



W 



W 



W 



q + hei 

q z q + 2h€i 
q + 

q + hei 
q z q + h(ei + ej) 
q + hei 

q + hei 
q z q + h(ei + 6j) 
q + hej 



(z + h) 
9(h) ' 



+ gtf) 

6(z) 9(h + gij ) 
0(h) e(qy) 



i + J) 



where q^ =< q,€i — €j >. 

The relation between R B and the face weights is given by 

^R B (z - ^f<P(z)t^,<PHV£ +lm) r = 



q + he k 
q z - w q + h(e k + e r 
q + ht s 



(4.6) 



(4.7) 
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In what follows we use the concise notation 



W?[k + m] = W 



q + he k 
q z -w q + h(e k + e m ) 
q + he s 



Then the dual relation to (|4.7Q is 

E 



,v)l +hlkj 'Mz)ltl { t + ~ em)i 'R B (z - w)%, = W^[k + m}^z)l +hl ^(w) q+H " k+ ~ em)j 



>q+hc k 



>q+he s 



(4i 



In |21| another L-operator L appeared. It is related to L in the following way: 



(4.9) 



i'j< 



fl(z + 7 + qjj) j-j 0(7 + q nj ) m^. 



9{z) 



0(q n i 



In what follows we need to remove from the quantum L-operator algebra ( 3.27|) the non- 
holomorphic dependence on the spectral and quantization parameters. This can be achived by 
considering the following transformation of the L-operator 



L{z) -> e a{z)Q e-P Q L{ Z y Q e- a{z)Q , 



(4.10) 



where a(z) is an arbitrary function of the spectral parameter and (3 is a complex number. Since 
the transformed L-operator also has the structure WP, then the following formula is valid 



L 2 (w)e a{z)Ql = e a(z)Ql L 2 (w)e ha(z)r °, 



(4.11) 



where the notation r = J2i En <S> En was used. 

Recalling that the L-operator ( |3.35| ) satisfies the quantum L-operator algebra ( |3.27| ) and 
using eq. (|4.11|) one can easily establish the algebra satisfied by the transformed L: 



Ri 2 (z,w)L 1 (z)R 21 (w)L 2 (w) = L 2 (w)Ri 2 (z)Li(z)R( 2 (z - w), 
where the matrices R, R and R F are 



R 12 (z,w) = e -^)Q^ w )Q*+PQ*R 



i2(z,w)e 



a{z)Q 1 +a(w)Q 2 -0Qi 



R 12 ( Z ) = e ha(z)r -a(z)Q 1+ /3Q 2Ri2 ^ e 



a(z)Q 1 -/3Q 1 



(4.12) 
(4.13) 

R( 2 (z,w)= (4.14) 

e -a(z)Q 1 -a(w)Q 2 +&{a{w)-a(z))r +PQi R F / w ^ e a(z)Q 1 +a(w)Q 2 + ^(a(w)-a(z))r -f3Q2 



Since the transformation in question keeps the form of the quantum L-operator algebra intact, 
the transformed matrices R, R and R F also satisfy all the compatibility conditions. In partic- 
ular, the transformation ( [4.14|) defines another solution of ( |3.25|) . For (3 = this was observed 
in M . 
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For the particular choice /3 = 2ni^—j an d ct(z) = liti^-^ + P we find that the matrices 
R, R and R F are given by the same formulae ( |3.13| ), (|3.16|) , ( |3.24j ) with the total change of 
$(z, s) = q^(s) ® '(0) e27rt7r? by the meromorphic function §^jfk ■ The transformation with 
such a choice of a and /? transforms (up to an unessential multiplier) the L-operator (|3.35|) into 

_ e{z + q %J + 1 ) n %nj + 7) 



which is a quasi-periodic meromorphic matrix function of the spectral parameter: 
L{z + 1) = L(z), L(z + t) = e- 2m ^ +n) e- 2mQ L(z)e 2niQ . 

he-2- 

We assume that the L-operator is of the form WP, where Pj = e ld n. The replacement of 
P by P preserves all the consistency conditions because the P-matrices depend only on the 
difference qi — qj. Thus, eq. (|3.27|) with P-matrices defined via &(z,s) = o^(s) Te ^ eTS to the 
meromorphic version of the quantum L-operator algebra while (|4.15|) provides its particular 
meromorphic representation. In what follows we use only this meromorphic version. 

Comparing ( 4.15| ) to ( [4.9| ) we can read off that L and L are related in the following way 

Since the combined transformation ( |4.9| ), ( [4.16 ) from L to L depends only on q we can conjecture 
that any representation L of the quantum L-operator algebra (|3.27f) is gauge equivalent to some 
representation L of (fO| ) with a gauge-equivalence defined as 



=^z)^(z)^ s ^^ L ffr (z). (4.17) 

Now we are in a position to prove this conjecture. Suppose L be an abstract L-operator 
satisfying algebra ( |4.1| ) and introduce L by eq. ( |4.9|) . Assume that L has the structure WP, 

where the entries of the diagonal matrix P are P j = e 1 s<3 * and the entries of W commute with 
qi. Then substituting L expressed via L in Q4.1D and performing the straightforward calculation 
with the use of (|4.8|) one finds an algebra satisfied by L: 



i'jl 

abed 



abed 



Performing the summation in %' and f with the help of (|4.4[) we obtain 

TT/6r , , , j. #(w + q ac + M a6 - Kdjc) -,-r %„ c + M nfe - h5 jc ) ~ 

g ^ + .] w. ^ n g^+Mu-^ ^w^w 

jjrj r . , 7l , 0(z + q ib + M ik - M bc ) „ 6(q nb + %8 nk - %5 bc ) ~ , . ~ , . 

E "2b + fl W ^ n o (qm + u vk -M, t) L ^ )L ^- 
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Let us introduce an operator L by inverting ( 4.16|) . Substituting this L in the last formula, 
taking into account the nonzero components of the face weights and multiplying both sides by 
the function 

Un^k Q[Qnk) Iln.^i #(<?m + ftS n k ~ ftSik) 
Un^j %nj) Tln^l %n* + ft^nj - HSji) 

we finally arrive to the algebra satisfied by L: 

cfc 6>(w + q ks + ft - M js ) Iln^fc 1 



s 
s 
s 



hS js ) 



0{w) 



L kj (z)L sl (w) + 



I 0(W + q is - ftS js ) Un^i 9(<lns + hb nk - ftS js ) 
Ojqki + ft) 0{W + q ks - M js ) Iln^fc + fi$ni 



L k j(z)L s i(w) + 



ft8 



0{w) Unjts (tins + ft$nj - ftSjs 

(z + q is + ftS ik - ftSis) 



ki — ft) 

' e( qij + ft) 
Aiij - *) 



Lij(z)L a i(w) 



26,, 



x 



Un^i 9(c[ns + ftSnk ~ ft$ls 



Un^s 9{<ln S + ftSnl ~ ft$ls 



-L kl (w)L sj {z) + 



0(z + q is + ftdik - ftSjs) Ylnjki 0(q ns + ft$nk - ftSjs 



-L kj {w)L sl (z). (4.18) 



6{z) Un^s 9{<lns + ftSnj ~ ftSjt 

Here in the second and the third lines i ^ k. The ratio of products of theta-functions occurring 
in each term in ( 4.18 ) allows one to take off the sum over s, e.g., when i ^ k, we have 

Tln^i %ns + ft^nk ~ ft^js) _ 



Lln^s ®{<lns + ft$nj ~ ft$js 
7 %fci) , 



Si, s. 



0{q ki + ft)B{q 



jo 



i.i 



0(q ki - ft) 



m + 



8 k ,(l — 8, 



kjj 



Si 



6(h) 



+ 



)6(h) 



S, 



\ 13 0(q ik + ft) ' 6(q jk + ft)6(q ik/ 
e{q k j)6{ft) 



0(q kj - ft)6{qji + ft) 



To compare (ELT8D to ($27\) we rewrite relation (ggg) with the help of eq. flOiD in the 
following form: 



w)P 1 R 21 (w)Pr 1 L 1 {z)R 21 {w)L 2 (w) 



P 2 R 12 (z)P 2 1 L 2 (w)R 12 (z)L 1 (z)R( 2 (z 



w). 
(4.19) 

In the component form algebra ( |4.19| ) is presented in Appendix C. Comparing the components 
of ( |4.18| ) to the ones of ( 4.19|) we establish that they coincide up to the overall multiplicative 
factor 6(z — w)6(h) 2 . Thus, we have shown that any representation of algebra ( p.27|) by the 
transformation ( |4.17| ) turns into a representation of (ITT). The connection established gives 
right to assert that algebra ( 3.271) possesses a family of iV-commuting integrals and that the 
formula of the determinant type ( 3.4(J| ) for the commuting family proved in ]2T| is also valid for 
the L-operator ( |3.35|) . 
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5 Conclusion 



In this paper, we described the dynamical .R-matrix structure of the quantum elliptic RS model. 
The quantum L-operator algebra possesses a family of commuting operators. It turns out that 
this algebra has a surprisingly simple structure and can be analysed explicitly in the component 
form. Furthermore, one can hope that the problem of finding new representations of the algebra 
obtained is simpler than the corresponding problem for the algebra RLL = LLR. 
There are several interesting problems to be discussed. 

First, we recall that in the classical case we obtained two different Poisson algebras, which 
lead to the same classical L-operator algebra. Only one of them was quantized. It is desirable 
to quantize the second one and to show that the corresponding quantum L-operator algebra is 
isomorphic to the algebra obtained in the paper. 

The elliptic RS model we dealt with in the paper corresponds to the Atv-i ro °t system. 
It seems to be possible to extend our approach to other root systems and to derive the corre- 
sponding L-operator algebras. To this end, one should find a proper parameterization of the 
corresponding cotangent bundle. 

Generalizing our approach to the cotangent bundle over a centrally extended group of 
smooth mappings from a higher-genus Riemann surface into a Lie group, one may expect 
to obtain new integrable systems. 

It is known that the CM systems admit spin generalizations [34, 35, 3(J. Recently, the spin 



generalization was found for the elliptic RS model |37|. However, the Hamiltonian formulation 



for the model is not found yet. One may hope that in our approach the spin models can arise 
as higher representations of the L-operator algebra. 

Probably, the most interesting and complicated problem is to separate variables for the 
quantum elliptic RS model. Up to now only the three-particle case for the trigonometric RS 
model was solved explicitly |5S| . One could expect that the L-operator algebra obtained in the 
paper may shed light on the problem. 
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Appendix A 



In this Appendix, we prove that the limiting r-matrix ( p. 20 ) satisfies the classical Yang-Baxter 
equation. For this, let us write the equation, which follows from the Jacoby identity for the 
bracket (E18D : 



[rl v 2 {z, w), rl p 3 (z, s) + r%{w, s)} + [r^(z, s), r%(w, s)} 
+-T 3 P -\8){Tg(8),*%(z,v,)) + -IT -\z){Tt(z),v^(w,s)] 



(A.L 



w){T2(w),v%(z,s))=0. 



Here 
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+ pi®i Xs H?H?r${z,w) +pti s H?H?r%{z,w) 
+ p 2 tr u H^HTr £ 3 2(z,w)I®I. 



(A.2) 



To study (|A.1|) , one needs to know the Poisson bracket of T £ (z) with H vp . This bracket can 
be easily derived from eqs. (|2.13| ), Q2.14|) , and (|2.17p : 



\t\ -\z){Tl(z), Hf } = E H?H?*(e - p, g y )(i& - h) ® £ 
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(A.3) 



E Ht P Hf^z - p, qij )E l3 ® Eji + E HrHfHY${z - rj, Qij )E i: 

i¥=3 i+3 

E (He -p,Q) + $(v -e,0) + *(z - r,, 0) - S(z - p, 0)) (J5y - 1/) ® 



Eq. (|A.l ) holds at arbitrary values of all the parameters. Without loss of generality one can put 
p = and f] = ae, where a and e are real. Let us now perform the change of variables Hf 1 : 

H* v = l + h? 

and consider the expansion of the l.h.s. of (|A.l ) in powers of e and h. 

Let us note that the matrix r ev (z, w) has the following expansion in powers of h at e — rj — > 0: 



r e,? (z, io ) = r(z, w) + r r eg(z, to, e)h + r (2> (z, w, e)h 2 + o(e — rj) 



(2), 



(A.4) 



where r(z,w) is given by Q2.20Q . The matrix r r eg(z, w, e) is regular when e — > and the e- 
dependence of r^ 2 ^(z,w,e) is inessential. Substituting (|A.4j ) into the bracket T _1 {T, r} one 
gets: 

T- J {T, r} = rSg^T-^T, /»} + rC^T^T, /»}. (A.5) 

Clearly, eq. (|A.l ) should be satisfied in any order in h and e. Since we are interested in finding 
an equation for r, we consider only the terms independent of h and e in the expansion of 
eq. flA.l ). The terms of zero order in h and e occurring in ( |A.1| ) come from r Er) and from the 
first term in ( A.5|) . However, from the explicit expression for the bracket {T,H} one can see 
that 

{T,h}\ h=0 = o(e) 

Now, taking into account that r r eg(e) is regular at e — > 0, we conclude that the last three terms 
in ([A.l|) do not contribute. Thus, in the zero order in h and e, eq. (|A.l|) reduces to the CYBE 
for r(z, w). 

The remarkable thing is that the main elliptic identities (see Appendix B) follow from 
the Jacoby identity for the bracket ( |2.18 ) or, equivalently, from the Yang-Baxter equation for 
r(z, w). 
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Appendix B 

Here we present the basic elliptic function identities, formulated as a set of functional relations 
on $(z,w) @: 

x)$(w, y) = x - y)$(z + w,y) + $(z + w, x)$(w, y — x), (B.l) 
$0, x)$(z, y) = x + y) ($(z, 0) + 0) + $(y, 0) - $(z + x + y, 0)) , (B.2) 
$(z -w,a — b)$(z, x + b)$(w, y + a) - $(z —w,x — y)$(z } y + a)$(w, x + b) = 
<f>(z, x + a)<f>(w, y + b) ($(a - b, 0) + $(x + 6, 0) - $(z - y, 0) - $(a + y, 0)) . (B.3) 



Eq.( p.2|) is the limiting case of eq. (|B.l|) where w — ► 2, and eq. (p.3| ) is a consequence of ( B.l| ) 



and ( |B.2| ). Note that the exponent term e 2 ^( 1 / 2 ) zs+27r * s T-T i n $>(z, s) as well as the linear term 
— 2£(l/2),z + 27ri^5f m "^(^jO) are irrelevant since they drop out from ( B.1 - B.2]) . 



To establish the unitarity relation for R, one also needs the identity involving the Weierstrass 
V-f unction 

$0, s)Q(z, -s) = V{z) - V(s), 
and to prove eq.( [2.50|) , the following relation between the derivatives of $ is of use: 

d${z,qij) _ d<5>{z, qij ) 



dz dqi 



Appendix C 

In this Appendix we present the quantum L-operator algebra 

Rf 2 {z-w)P 1 R^{w)P^ 1 L 1 {z)R 21 {w)L 2 {w) = P 2 R^(z)p- 1 L 2 (w)R 12 (z)L 1 (z)R^(z-w) (C.l) 
in the component form. The l.h.s. of ( |C1[ ) has the form 

Qki)®{h, Qik)®(h, Qkj ~ h)L ij (z)L k i(w)E ij <g> E u 

i^k;j,l 

- $ (^> Qu)®(h, Qij ~ h)$(w, q kj - h)Lij(z)Lji(w)Eij <g> E M 

i^k;j,l 

+ J2 $ (^' Qki)&(w, Qki)®(h, % ~ h)L ij (z)L a (w)E ij ® E M 

i^k;j,l 

+ ^ $(z- w, q ik )®(h, Qij - h)L kj (z)L il (w)E ij ® E H 

i^k;j,l 

^ $(2: - w, qik)${h, qkj - ti)$(w, qij - h)L kj (z)L j i(w)E ij ® £ fci 
+ 51 ®(z - w,q ik )$(w,q ik )$(h,q kj - h)L kj (z)L k i(w)Eij ® E kl 

i^k;j,l 

+ ®( z ~ w > h)<$>(h, q kj - h)L kj (z)L M (w)E kj <g> E kl 

j,k,l 

$0 - w, h)$(w, h)$(w + h, q kj - fi)L kj (z)L jl (w)E kj ® E kl 

j,k,l 
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The r.h.s. of ( |C.1|) reads 



J2 $0, Qki)H% Qu ~ ti)$(h, qi^Lki^Lij^Eij ® E kl 
+ *(^' Qki)$(h, Qij ~ h)$(z - w, q tj + h5i j )L kj (w)L u (z)E ij ® E kl 

i^k;j,l 

Qki - K)&(z, Qu - qi^Lk^Lij^Eij ® E M 

Qkj - h)$(z, q-ij - K)$(z - w, qij + h5i j )L kj (w)L jl (z)E ij <g> E H 

i^k;j,l 

+ 9ik)&(h, qu - h)$(h, q lj )L kl (w)L kj (z)E ij <g> E w 

+ ^ q ik )®(h, q kj - h)$(z - w, % + h5i j )L kj (w)L kl (z)E ij ® £7 W 
+ $ ( z > fti - qi j )L u (w)L ij (z)E ij <g> £;« 

+ J! 9y _ _ ^' ®i + ^i j )L ij (w)L il (z)E ij ® 

5] $(z, ft)$(z + ft, <fe - qi^LuiwjL^Eij <g> E a 

51 ^) $ ( 2 + ^ ftj ~ ~ w > *i + fiS j i)L ij (w)L j i(z)E ij ® £^ 
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